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Abstract.  We  analyze  strain  induced  effects  on  the  capacitance- voltage  (C-V)  profile  for  MBE 
grown  GaAs/Ino.iGao.9As/GaAs  HFET  structure.  The  calculations  of  C-V  profile  were  made 
using  a small-signal  approach  developed  for  the  quantum  well  (QW)  structures.  The  self-consistent 
numerical  simulations  and  results  of  measurements  show  that  strain  causes  significant  changes  in 
the  electron  distribution  of  the  QW  structure  and  its  C-V  profile. 


1.  Introduction 

Capacitance-voltage  (C-V)  profiling  is  on  of  the  methods  for  the  investigation  and  design 
of  heterostructures  quantum-well  (QW)  electronic  and  optoelectronic  devices  [1-5],  This 
technique  is  used  for  determining  band  offsets  at  the  heterointerface.  In  the  case  of  the 
psedomorphic  structures,  in  particular,  the  GaAs/Ino.iGao.9As/GaAs  system,  with  suffi- 
ciently thin  central  layer  the  lattice  constant  mismatch  is  accommodated  by  internal  strains 
rather  than  by  formation  of  misfit  dislocations.  The  strains  lead  to  two  consequences:  the 
shifts  of  the  conduction  and  valence  bands  in  the  well  for  the  structures  grown  in  [100] 
direction  [6]  and  the  built-in  piezoelectric  field  for  the  structures  grown  in  [11 1]  direction 
[7].  However,  with  the  exception  of  the  Ref.  [2],  the  calculations  of  band  offsets  from 
C-V  experimental  data  were  made  without  the  strain  effect.  At  the  same  time,  the  strain  is 
sufficient  large  and  produces  the  important  changes  in  interband  optical  transitions  [6].  The 
purpose  of  this  work  is  to  study  the  strain  effect  in  GaAs/Ino.iGao.9As/GaAs  quantum- 
well  HFET  structure  using  the  self-consistent  numerical  simulations  and  to  compare  the 
theoretical  results  with  the  measured  C-V  profile.  The  new  simulation  technique  of  the 
C-V  characteristics  for  QW  stmctures  based  on  the  small-signal  approach  [8]  is  developed. 

2.  Experiment 

Experimental  sample  was  grown  by  molecular  beam  epitaxy  using  a Riber32P  system. 
Undoped  GaAs  layer  with  thickness  0.5  mm  was  used  as  a buffer.  The  conductive  channel 
consists  of  n-GaAs  (10  nm)  followed  by  n-InGaAs  (25  nm)  and  n-GaAs  (75  nm)  layers. 
Structure  was  capped  by  n-GaAs  contact  layer.  The  channel  composition  was  optimized  to 
obtain  high  linearity  of  the  static  and  dynamic  HFET  characteristics.  The  process  sequence 
for  device  fabrication  is  listed  bellow: 

1.  Source  and  drain  AuGa/Ni/Au  contacts  formation; 

2.  Device  isolation  by  mesa  etching  and  following  proton  implantation; 

3.  Gate  electrode  formation  (E-beam  lithography,  gate  recess  and  Ti/Pt/Au  metallization); 

4.  Device  passivation  and  contact  pads  metallization. 
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C-V  profile  for  Schottky  diode  test  structures  ( 1 00  mm  in  diameter)  was  obtained  from 
capacitance-voltage  measurements  at  1 MHz. 

3.  Theoretical  model 


The  C-V  (or  apparent)  profile  Nap  (W)  is  given  as  usual  by: 


Nap  = 


C\V) 
ee(dC/dV ) 


, W(V)  = 


e 

COO  ’ 


(1) 


where  e is  the  electron  charge,  e is  the  semiconductor  permittivity,  and  V is  the  reverse 
bias  voltage  applied  to  the  Schottky  barrier.  Capacitance  C(V)  can  be  calculated  as  [8]: 


C(V)  = 
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where  z is  the  space  coordinate,  L is  the  total  structure  thickness,  and  A n(z,  V)  is  the 
perturbation  of  electron  concentration  due  to  the  small  voltage  variation  A V,  p is  the  dis- 
tribution function  for  the  perturbed  space  charge  density.  Usually  the  perturbation  A n(z,  V) 
is  calculated  as  the  difference  beetwen  two  static  states  «(z,  V + A V)  and  n(z.  V)  [3,  4]. 
Such  approach  demands  the  high  accuracy  and,  consequently,  the  considerable  computer 
resources.  Recently  we  reported  about  very  effective  simulation  technique  of  the  barrier 
capacitance  for  the  classical  structures  based  on  the  small-signal  approximation  [8],  Now 
we  applied  this  technique  to  the  QW  structure.  In  this  case  the  electron  concentration 
distribution  n(z,  V)  must  be  determined  as  the  sum  of  the  3D  and  2D  concentrations,  this 
later  is  obtained  by  following  expression: 


m 

n2D(z,  V)  = ^|4Mz,  V)|2»/(V),  (3) 

i=l 


where  m is  number  of  bound  states  in  QW,  /?,-  is  the  electron  occupation  of  state  with 
eigenenergy  Et  and  wavefunction  The  static  2D  electron  distribution  is  determined 
from  the  Schrodinger’s  equation.  The  Hamiltonian  H in  the  effective-mass  approximation 
for  the  conduction  band  can  be  written  as  [2]: 


H2  d a2(z)  d 1 
2 a(z)  dz  m(z)  dz  a(z)  L' 


(4) 


were  a(z ) and  m (z)  are  the  position-dependent  lattice  constant  and  effective  mass.  The 
effective  potential  energy  Vef  includes  four  terms: 


Vef(z) 


Ec(z ) + 2a  c 


+ eU(z)  + Vxc(z), 


(5) 


were  Ec(z)  is  the  conduction-band-edge  energy,  ac  is  the  hydrostatic  deformation  potential, 
Cn  and  Ci 2 are  the  elastic  stiffness  constants,  sxx  = syy  is  the  strain  in  the  plane  of  the 
epitaxial  growth  [001],  Vxc(z)  is  the  exchange-correlation  energy.  The  Hartree  electrostatic 
potential  U (z)  is  determined  using  the  one-dimensional  Boltzmann-Poisson  equation.  So, 
the  calculation  of  the  static  electron  distribution  includes  the  quantum  effects.  After  this, 
the  barrier  capacitance  and  apparent  profile  are  simulated  in  the  classical  small-signal 
approximation  [8].  The  comparison  of  such  “semi-classical”  technique  with  the  usual 
quasistatic  approach  gives  that  capacitance  error  is  bellow  0.2%. 
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4.  Results  and  discussion 

In  Fig.  1 are  presented  the  simulated  3D  and  3D+2D  electron  distributions  at  the  reverse 
bias  V — 0 V.  One  cane  see  that  the  quantum  effect  changes  the  electron  profile  in  the  QW 
ofHFET. 
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Fig.  1.  Conduction-band-edge  energy,  lowest  three  sub-bands  and  electron  distributions  in  QW. 

Next  Fig.  2 shows  the  results  of  C-V  profiling  and  simulated  apparent  profiles  for  three 
cases:  (i)  unstrained  QW;  (ii)  strained  QW,  the  static  electron  distributions  are  calculated 
in  classical  approach;  (iii)  strained  QW,  the  self-consistent  simulation  of  the  static  electron 
distributions.  On  can  conclude  that  strain  effect  makes  the  considerable  influence  to  the 
apparent  profile  by  the  shift  of  conduction-band-edge  energy  in  the  QW.  At  the  same 
time,  the  inclusion  of  the  size  effects  produces  the  lesser  influence  for  the  given  structure. 
The  some  discrepancy  between  the  right  slope  of  the  simulated  apparent  profile  and  the 
experimental  data  cane  be  explicated  by  the  effects  of  the  edge  capacitance  and  increasing 
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Fig.  2.  Experimental  and  simulated  apparent  profiles. 
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Fig.  3.  Extracted  mole  fraction  and  doping  profiles. 


of  the  channel  resistance.  The  extracted  mole  fraction  and  doping  profiles  are  presented  in 
Fig.  3. 

The  obtained  thicknesses  of  the  first  GaAs  n-layer  is  48.5  nm  and  the  second  n-layer 
is  9 nm,  the  thickness  of  the  QW  is  24  nm.  The  doping  level  of  the  GaAs  n-layers  is 
5.6  • 1017  cm-3.  This  results  demonstrate  the  validity  of  the  proposed  “semi-classical” 
small-signal  technique  for  the  C-V  profiling  simulation  of  the  real  Q W structures  including 
strain  effects. 
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